Frequency and Q-factor control of nanomechanical resonators 
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We present an integrated scheme for dielectric drive and read-out of high-Q nanomechanical 
resonators which enables tuning of both the resonance frequency and quality factor with an applied 
DC voltage. A simple model for altering these quantities is derived, incorporating the resonator's 
complex electric polarizability and position in an inhomogeneous electric field, which agrees very well 
with the experimental findings as well as FEM simulations. By comparing two sample geometries 
we are able to show that careful electrode design can determine the direction of frequency tuning 
of flexural in- and out-of-plane modes of a string resonator. Furthermore we demonstrate that the 
mechanical quality factor can be voltage reduced more than fivefold. 



Control of small-scale mechanical systems is essen- 
tial for their application. Resonant micro- and nano- 
electromechanical systems (M/NEMS) have both proven 
themselves technologically viable (frequency filtering in 
cell phones [T], gyroscopes [2J, atomic force microscope 
(AFM) cantilevers[5]) as well as shown great promise for 
next-generation sensor applications (mass sensors [IHS], 
resonant bio sensors [7] and ultra sensitive force sensors 
[H105])- Three areas of development are central to real- 
izing the potential of high performance resonant micro- 
and nanomechanics: advancement of high Q geometries 
and materials; improved readout schemes for mechani- 
cal motion, including compactness and integr ability; and 
increased control of the resonant behavior of the me- 
chanics. In the field of nanomechanics, the last years 
have seen the advent of high Q silicon nitride strings un- 
der high tensile stress [TU], as well as indications that 
the stress rather than the nitride is responsible for the 
high mechanical quality factor [11) . Subsequently, effi- 
cient integrated drive and read-out schemes have been 
developed to detect the sub-nanoscale motion of small- 
scale resonant mechanics [T2]. Very good tunability of 
the resonance frequency can be achieved by capacitive 
coupling of the nanomechanical element to a side elec- 
trode [13 . However, the required metalization of the 
resonant structure reduces the room temperature qual- 
ity factor significantly [3j via Ohmic losses. In our lab, 
an efficient, room-temperature microwave mixing scheme 
has been developed for readout [15] . as well as a dielec- 
tric drive mechanism to actuate mechanics regardless of 
their material make-up[16j. importantly obviating the ne- 
cessity to metallize otherwise low-loss dielectrics. In this 
paper we present a continuation of this development that 
enables tuning of both the frequency and quality factor of 
nanomechanical resonators in the context of this highly 
applicable and integrable scheme. Using the combined 
dielectric actuation and microwave readout schemes, we 
theoretically develop the means to controllably raise and 
lower the resonant frequency of various flexural modes 
of our mechanics, as well as to broaden the mechani- 
cal resonance linewidth. The latter so-called Q control 
[T7][TH]is widely used in AFM measurements and grants 



faster image acquisition rates by decreasing the mechani- 
cal response time. In the presented scheme this Q control 
is achieved by simply applying a DC voltage instead of 
the normally employed active feedback mechanism and 
is thus well suited for integration. The theoretical re- 
lationship between the design of the electrodes and the 
resulting control of a given mode is validated both by 
experiment and simulation. 

Our system is depicted in Fig. [T] A nanomechanical sil- 
icon nitride string is situated between a pair of near-lying 
electrodes (Fig.JIJa)). A DC voltage across the electrodes 
produces an electric field that induces a dipolar moment 
in the silicon nitride string, which thus experiences a 
force in the gradient field. This force can be modulated 
by an additional AC voltage, actuating the mechanical 
resonance |16|. Additionally, the mechanical resonator is 
dielectrically coupled to an external microwave cavity via 
the same electrodes. An equivalent circuit diagram is 
shown in Fig. [TJb) . Deflection of the beam translates 
into a change of the capacitance C m (t) between the two 
electrodes and thereby modulates the cavity transmission 
signal. The mechanical oscillation can then be detected 
by demodulating this signal. A more detailed description 
of this room temperature heterodyne detection includ- 
ing the demonstration of self-oscillation caused by cavity 
backaction can be found elsewhere [T5] , 

To enable direct actuation of the mechanical resonator, 
we introduce a microwave bypass between ground and 
one of the electrodes using the single layer capacitor 
(SLC) 19J Cb y ■ Thus a DC and RF voltage can be applied 
to this electrode, whereas the other electrode is grounded 
via the microstrip cavity (compare Fig. [ijb)). Conse- 
quently, the presented integrated configuration combines 
the above mentioned cavity-enhanced detection [15] with 
dielectric driving and tuning [16]. Note that the drive 
and detection can be realized with a single set of elec- 
trodes, since crosstalk is avoided by operating the cav- 
ity at 3.5 GHz, well-separated in frequency from the me- 
chanical fundamental flexural mode frequency of here 
6.5 MHz. 

For this study, two sample geometries for obtaining op- 
timized gradient field coupling are fabricated from high- 
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FIG. 1. (color online), a) False color SEM micrograph of a 
55 um long silicon nitride resonator in the configuration de- 
picted in c). The gold electrodes forming the capacitor C m (t) 
are highlighted in yellow, the silicon nitride resonator and its 
clamping points in green, b) Equivalent circuit diagram of the 
transduction scheme with an inductively-coupled microwave 
cavity - represented by the capacitance C and inductance L - 
for dielectric readout. C m (t) is the capacitance between the 
gold electrodes which is modulated by resonator displacement. 
The microwave bypass capacitor Cb y allows the additional ap- 
plication of a DC and RF voltage, U and Urf respectively, and 
thus dielectric tuning and driving. c,d) Schematic cross sec- 
tion with simulated field lines for the elevated and lowered 
geometry. The cross section of the silicon nitride beam is 
260 nm x 100 nm. The arrows in c) describe the directions of 
the in-plane and out-of-plane oscillation. 



stress silicon nitride using standard electron beam lithog- 
raphy, reactive ion etching and subsequent hydrofluoric 
wet etching to obtain suspended strings. To minimize 
Ohmic losses of the microwave circuit originating in the 
bulk chip supporting the mechanical resonator, we use a 
fused silica substrate instead of silicon. The geometries 
are schematically shown in Fig. [ljc,d). Referring to the 
beam's position with respect to the electrodes, the two 
structures will from now on be referenced as "elevated" 
(Fig.fiTc)) and "lowered" (Fig. HT d)). The centerpiece of 
each structure is the 55 um long silicon nitride string res- 
onator with a rectangular cross section of width 260 nm 
and height 100 nm. The freely suspended resonator is 
bordered by two vertically offset gold electrodes, one of 
which is connected to the microstrip cavity with a res- 
onance frequency of 3.5 GHz and a quality factor of 70, 
while the other electrode leads to the SLC. The essential 
difference between the geometries is the vertical position- 
ing of the beam with respect to the gold electrodes. This 
affects the dielectric environment and thereby the electric 



field lines as depicted in Fig. []Jc,d). In the elevated de- 
sign shown in Fig.[ljc), the upper edge of the electrodes 
coincides with the lower edge of the beam, whereas in the 
lowered design in Fig. |TJd) , the upper edge of the beam 
coincides with the lower edge of the electrodes. The simu- 
lated electric field lines for both geometries are obtained 
from finite element simulations using COMSOL Multi- 
physics and allow us to extract the electric field along the 
x- and y-direction. These inhomogeneous electric fields 
cause force gradients for the in- and out-of-plane modes 
of the resonator. They thus alter the restoring force of 
the respective mode and thereby its resonance frequency 
|16) . At the same time, the mechanical quality factor can 
be altered with the DC voltage, as the strong electric 
field and high field gradient lead to velocity-dependent 
dielectric losses in the beam material. This frequency 
and lincwidth tuning can be described by a simple model 
which agrees very well with our experimental findings 
and finite element simulations. The resonance frequency 
can be tuned over 5 % and the resonance linewidth can 
be increased by a factor of 6. 

We find the force gradient to be proportional to the 
square of the voltage and thus expect a quadratic depen- 
dence of the resonator resonance frequency on the ap- 
plied DC voltage. This can be derived from the energy 
of the induced dipolar moment p — aE of the dielectric 
resonator in an external electric field E, where a is the 
polarizability of the silicon nitride beam. Assuming that 
p and E are parallel, i. e. using a scalar, complex polar- 
izability a = a' + ia" and introducing a dependence of 
the electric field on the variable coordinate £, the energy 
W reads 

W = p ■ E = pE = aE 2 = (E(£)f(a' + ia") (1) 

Here £ can be the x- or y-coordinate (compare 
Fig. [ijc)), so the following considerations apply to both 
the in- and out-of-plane mode. Assuming E(£) = Eq + 
Ei£ for small displacements, the total energy can be sep- 
arated into a real (stored) and an imaginary (dissipative) 
part. 

Wstored = «' {El + 2E E^ + Eli 2 ) (2) 

W loss = a" (E 2 + 2E E 1 £ + E 2 £ 2 ) (3) 

The stored energy Motored is related to the frequency 
shift of the mechanical resonator, as the second derivative 
of this energy provides an additional force gradient i. e. 
an electrically induced spring constant k e : 

8F e d Wgtored /„2 r A \ 

= = uE 1 (4) 
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Hence the shift in resonance frequency caused by the elec- 
trical spring constant k e can be expressed in experimen- 
tally accessible units as 

f3 2 U 2 a' 
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DC tuning voltage U-U (V) DC tuning voltage U-U (V) 

FIG. 2. (color online). Quadratic tuning of the mechanical resonance frequency with DC voltage U for the two different 
geometries. The graphs show the deviation of the resonance frequency / from the natural resonance frequency fo of the 
resonator's respective mode (in- or out-of-plane, depicted as open and filled symbols respectively) for different microwave 
cavity pump powers indicated in dBm. The solid lines are a fit of the model, (a) The force gradient has the same parity for 
the in-plane- (open symbols, small curvature) as well as the out-of-plane mode (filled symbols, large curvature), (b) Tuning 
direction differs for the two modes. The out-of-plane mode tunes upwards, the in-plane-mode downwards in frequency. 



which assumes both, the electrical frequency shift 
k e /2mfo to be much smaller than the mechanical fre- 
quency /o and a geometry-dependent proportionality be- 
tween the applied voltage and the resulting electric field 
gradient E\ — f3U. Moreover, as predicted by our finite- 
element simulations the sign of the gradient (correspond- 
ing to the sign of /?) depends on the chosen geometry such 
that the out-of-plane mode changes its tuning direction 
between the elevated and the lowered design, which does 
not occur for the in-plane modes. 

The quadratic tuning behavior with DC voltage is 
found to agree very well with the experimental data, as 
displayed in Fig. [2] All measurements are conducted at 
room temperature and a pressure of 10~ 4 mbar. For each 
mode and geometry the mechanical spectrum of both 
modes is taken for different DC voltages and microwave 
powers. The driving voltage U T { is kept constant in every 
measurement. The values for U T { lie within 80 uV and 
1 mV depending on the particular mode and geometry. 

A Lorentzian fit to each mechanical spectrum yields 
the resonance frequency and the quality factor of the res- 
onance for each parameter set. The resonance frequen- 
cies lie around 6.5 MHz and the highest quality factor is 
340 000 for the out-of-plane mode in the elevated design. 
Note that the tuning with microwave power is a result 
of the effective microwave voltage [T5], and so is analo- 
gous to the tuning with a DC voltage. Subsequently, we 
fit / = /o + c(U — Uo) 2 + c mw /7^ w to the tuning curves 
shown in Fig. [5] using the natural resonance frequency 
fo and two tuning parameters for the DC voltage U and 
the effective microwave voltage [/ mw , as the DC and high 
frequency polarizability might differ. We also introduce 
the DC offset Uq to account for a shift (typically less 



than 1 V) of the vertex of the tuning parabola away from 
zero DC voltage, which is most likely caused by trapped 
charges in the dielectric resonator material. As the in- 
fluence of the microwave field on static dipoles averages 
out, there is no such shift resulting from the microwave 
voltage U mw . Consequently, we can extract the tuning 
parameters for each geometry and oscillation direction. 
With increasing voltage U and for the elevated geome- 
try depicted in Fig. [5Ja), both the in- and out-of-plane 
mode tune to lower frequencies, whereas for the lowered 
design (Fig. [2jb) ) , the out-of-plane mode tunes to higher 
frequencies, as predicted by our simulations. The solid 
black lines in Fig. [2] show the fit of our model with a 
single set of parameters for each mode in excellent agree- 
ment with the data. In the case of opposite frequency 
tuning, the initial frequency difference of the in- and 
out-of-plane modes can be evened-out, which leads to an 
avoided crossing caused by a coupling between the modes 
|20j . As the data points in this coupling region deviate 
from normal tuning behavior, they have been omitted in 
Fig. [lb). 



Altering the DC or effective microwave voltage does 
not only shift the resonance frequency, but also influences 
the damping T — 2nf/Q = 2nAf of the mechanical res- 
onance and thereby the measured linewidth A/. Like 
the resonance frequency, the damping T varies quadrat- 
ically with increasing voltage. This can be understood 
by analyzing the dissipated energy Wi oss given by equa- 
tion A time average of this quantity over one period 
of mechanical vibration £(t) = £o cos(wt) gives 
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W^T S = — / [2E E^ Q cos(uit) + E^l cos 2 (c^)] dt 
o 

= \*"Eie (6) 

Here we omit the Eq term (as a"(cu = 0) = 0, other- 
wise static electric fields would lead to dissipation). As 
the mechanical stored energy Wmoch = ^towq£§ is much 
larger than the electrical energy Whored > one can approx- 
imate the additional electrical damping to be 



The ratio 



— = tan(0) = ^ (10) 

is then independent of all resonator parameters and 
can be determined from the two curvatures. 

The measured values for damping and tuning cur- 
vatures are cr = 5.2^ and c = 438^1, leading to 
tan(</>) = 0.037. By using the Clausius-Mossotti-Relation 
to first calculate the (lossless) a using e = 7.5, one can 
determine the dielectric loss tangent 



r e (£/) = 



& = P 2 U 2 a" 
2nW mcch 2-KmuJo 
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FIG. 3. (color online). Damping constant versus DC voltage 
for two different microwave powers, exhibiting a quadratic 
behavior. The solid lines are a fit of the model mentioned in 
the text. 

The measured damping versus DC voltage is shown in 
Fig. [3j It displays the quadratic behavior of the damping 
constant T = r + T e (U) = T + crU 2 of the out-of-plane 
mode in the elevated design for two different microwave 
powers. Here Tq is the intrinsic damping of the resonator 
and r e ([/) is given by equation Q. Again, the vertical 
offset between the two curves is explained by the effective 
microwave voltage acting analogously to a DC voltage. 
The solid lines in Fig. [3] are a fit of the model to the 
data, from which the curvature cr can be extracted. 

Using this curvature and equation Q, the imaginary 
part of the polarizability can be expressed by 

a = J2 ( 8 ) 

Similarly, employing the curvature c of the parabolic fre- 
quency shift and using equation ^ an expression for the 
real part a' can be given: 



tan((5) = — = 0.016 
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This value is well within the range of loss tangents 
reported for silicon nitride thin films [21] • Note that the 
time- varying capacitance C m (t) induces a dissipative cur- 
rent in the electrodes, which also leads to a quadratically 
increasing damping 113,. However, using values obtained 
from FEM simulations for the electrode capacitance and 
its variation with beam deflection [TO], we estimate that 
this damping is three orders of magnitude smaller than 
that caused by dielectric losses. The relevant effect for 
the additional damping with increasing DC voltage is 
thus the dissipative reorientation of the dipoles in the 
resonator caused by its motion in a static, inhomogenous 
electric field, as described by equation §6§. 

The DC voltage dependence of the mechanical damp- 
ing r was also measured at zero microwave power us- 
ing an optical detection technique [35] . The resulting T 
was within a few percent of the value extracted from the 
9 dBm curve in Fig. [3] demonstrating that a measure- 
ment at low microwave powers induces only negligible 
additional damping to the mechanical resonator. 

In conclusion, we show dielectric frequency tuning of 
over 5 % of the natural resonance frequency for nanomc- 
chanical resonators in an all-integrated setup that re- 
quires no metallization of the resonant mechanical struc- 
ture itself. This scheme thus maintains an excellent qual- 
ity factor of up to 340 000 at 6.5 MHz and 300 K. Fur- 
thermore, by careful design of the geometry, one can 
choose the tuning behavior of the out-of-plane mode to be 
either upward or downward in frequency and thus tune 
the two orthogonal resonator modes both in the same 
and in opposite directions. We demonstrate that dielec- 
tric losses become highly relevant when using nanoscale 
electrode geometries generating large field gradients pro- 
viding high tunability. This could be used as a Q control 
[T71 ITS! |2"3"] that does not require any active electron- 
ics such as a phase-locked loop but rather a single DC 
voltage. Such a Q control can be employed to increase 
the bandwidth of NEMS sensors significantly, leading to 
much more adaptable devices. Without the need for ac- 
tive electronics this could prove to be very well suited 
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for integrated designs. Together with microwave cavity 
backaction |15j the mechanical resonance linewidth can 
be controlled from a few Hz up to more than 100 Hz, 
thereby tuning the mechanical bandwidth by about two 
orders of magnitude. Finally, we imagine that the scheme 
presented can also be employed to build self-sensing AFM 
cantilevers [213 with tunable bandwidth and resonance 
frequency that are not subject to the bandwidth limi- 
tations of the normally employed piezo drive and could 
thus be used in multifrequency force microscopy schemes 
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